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Abstract. The effect of the four-spin cyclic exchange interaction at each plaquette 
in the S = 1/2 two-leg spin ladder is investigated at T — 0, especially focusing 
on the field-induced gap. The strong rung coupling approximation suggests that it 
yields a plateau at half of the saturation moment (m = 1/2) in the magnetization 
curve, which corresponds to a field-induced spin gap with a spontaneous breaking of 
the translational symmetry. A precise phase diagram at m = 1/2 is also presented 
based on the level spectroscopy analysis of the numerical data obtained by Lanczos 
method. The boundary between the gapless and plateau phases is confirmed to be of 
the Kosterlitz-Thouless (KT) universality class. 
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1. Introduction 

The multiple-spin exchange interaction has attracted a lot of interest in the condensed 
matter physics. Such many-body exchange interactions are realized in the two- 
dimensional (2D) solid 3 He [g§ and the 2D Wi erier solid of electrons formed in a silicon 
inversion layer ||, as well as the bcc 3 He [|J. Recently the four-spin cyclic exchange 
interaction, called ring exchange, has been revealed to be important even in the strongly 
correlated electron systems like the high-temperature cuprate superconductors. The 
analysis on the low-lying excitation spectrum of the d-p model indicated that the ring 
exchange should be taken into account in the simplified spin Hamiltonian describing 
the Cu0 2 plane 0. The study based on the Heisenberg Hamiltonian also revealed an 
evidence of the ring exchange in the Raman scattering spectrum . In fact such a four- 
spin interaction was derived from the fourth-order perturbation expansion of the square 
lattice Hubbard Hamiltonian with respect to t/U near half-filling [[/]. The neutron 
scattering experiment also suggested that the effect of the ring exchange appeared in 
the spin wave excitation spectrum of La2Cu04 ||. 

Recently the ring exchange has been supposed to be important also in the spin 
ladder systems. The significant difference between observed leg and rung bi-linear 
exchange coupling constants (Ji eg ~ 2J rung ) of Sr 14 Cu240 4 i was explained, assuming 
the existence of the ring exchange with the amplitude of about 14% of J rung J|. A four- 
spin exchange interaction described by a product of two-spin exchanges in a spin ladder 
was investigated by a field theoretical approach |L0| , where the possibility of a different 
type of massive phase from the Haldane phase is indicated [|TT|,|T2]]. in the nonmagnetic 
ground state. The recent density matrix renormalization group study suggested that 
increasing ring exchange constant J 4 brings about a quantum phase transition about 
J 4 ~ 0.3J mng for Ji eg = J rung However, the feature of the large- J 4 phase is still an 



open problem. We note that there exist integrable spin ladder models with four-spin 



exchange interactions [Tj,|T^]. But these models are somewhat differnt from our model. 




Figure 1. Illustration of a two-leg spin ladder with four-spin exchange interaction at 
every plaquette. 



As another interesting phenomenon caused by the ring exchange, a field-induced 
spin gap which should be observed as a plateau in the magnetization curve, was 
proposed in an S = 1/2 antiferromagnet on the triangular lattice JT7] . It was verified 
by the exact diagonalization ||18|| . In the spin ladder system, besides the original spin 
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gap []l9-2"Tfl, the magnetization plateau was also revealed to appear in the presence of an 
additional leg ||22|| , a bond alternation |23|,[24[] or some frustrated interactions p5|-[3T 



In the previous work |32| by two of the present authors (TS and YH), such a field 
induced gap due to the ring exchange in the spin ladder was investigated. In that work, 
the phenomenological renormalization study []33| suggested that a plateau appeared 
at half of the saturation magnetization in the spin ladder for a realistic parameter 
J 4 > (0.05 ± 0.04)J rung in the case of J\ cg = J mng . However, the critical point J 4c 
estimated by the one-magnon analysis did not agree well with that by the two-magnon 
one. The critical index at J 4c was also different from that of the Kosterlitz-Thouless 
transition ]34| predicted by the behavior of the energy gap. Thus some better studies 
are necessary to determine the critical point J 4c and the universality class for the 
transition between the plateau and gapless states due to the ring exchange in the 
spin ladder. In the present paper, the perturbation expansion from the strong rung 
coupling limit is performed to clarify the mechanism of the plateau formation and the 
quantum critical behavior around the phase boundary. Also, a precise phase diagram 
is presented, using the size scaling analysis based on the conformal field theory |35|-p7 



and the recently-developed level spectroscopy method [pq-|41 
also more generalized than that of the previous work. 



The parameter space is 



2. Spin Ladder with Ring Exchange 

We consider the S = 1/2 uniform antiferromagnetic spin ladder with the four-spin cyclic 
exchange at every plaquette, as shown in figure 1. The Hamiltonian of our model is 

H = Ho + Hz, 



H = J2(Ji + JiS 2 + J± 



2,3 J 



+ j 4 E(^ + p 4 - 1 ), 



(1) 



3=1 



H z = —H 22{S{j + S^j 



where J± and J± are the bi-linear leg and rung exchange constants, respectively. We set 
J_i_ = 1 throughout this paper. P 4 j is the cyclic permutation operator which exchanges 
the four spins around the j-th plaquette as Sij — > <Si,i+i > ^2,j+i "~ * Sij, 



P., 



PI 




(2) 



d a\ 
c &/' 

where the (1, 1) component denotes the spin state of Si^ and the (2, 2) component that 
of S2J+1- The strength of the four-spin ring exchange J 4 is assumed to be positive, as 
it is in the Cu oxides. The applied magnetic field is denoted by H . The magnetization 
of the bulk system is defined as m = M/L, where M = (J2j(St,j + ^fj))- We focus on 
a possible plateau at half of the saturation (m = 1/2) in the magnetization curve at 
T = 0. 



Field-Induced gap due to four-spin exchange in a spin ladder 



4 



3. Strong rung coupling approximation 



The condition of the quantization of the magnetization (16| predicted that Q(S — 
m) =integer is necessary for the appearance of the plateau at m, where Q is the period 
of the ground state and S is the total spin per unit cell. Thus the plateau at m = 1/2 
should require a spontaneous breakdown of the translational symmetry which results in 
the period of the two unit cells like the Neel state. 

In order to investigate the possibility and the mechanism of the plateau at m = 1/2 
in the spin ladder with the ring exchange, we consider the strong rung coupling limit 
(Jj_ 3> Ji, J4) |27|,iH a t first. The system is separated into isolated rung dimers when 
J\ = J a = 0. In this case, at m — 1/2, half of the rung pairs are triplet with S z = 1, 
and the remaining half are singlet. Since other possible selections have higher energies, 
we may restrict ourselves to these two states for each rung dimer and represent them 
by states | fr) and | Jj) of pseudo-spin S 




V2 ( 




(3) 



When J\ = J 4 = 0, each pseudospin S can take either | ff) or | JJ.) completely freely, 
which means high degeneracy of the m = 1/2 ground states. The introduction of weak 
J\ and J 4 resolves the degeneracy. The effect of a weak inter-dimer coupling can be 
taken into account by the perturbation expansion with respect to J\ and J4, resulting 
in the following matrix elements by straightforward calculation 



V 



I Mm) 

-(1/4)J! + 2J 4 




By comparing equation ([4j 
magnetic field H 

I Mj+i) 
f (1/4) J z -H + c 








-(1/4) ^ - J 4 
(1/2) Ji + J 4 






(1/2) J x + J4 
-(1/4) Ji - J 4 




I JJ-j-U-j+i) 
\ 



(1/4) Jx + J J 



(4) 



with an XXZ-type Hamiltonian including an effective 

-1) 



(jj-jjj-j+i 



V 





-{l/A)J z + a 
(1/2) 






(1/2) 3*» 

-(l/A)J z + a 




\ 




(l/A)J z + H + aJ 



(5) 



(6) 



we obtain the effective Hamiltonian for S 

n = E {J xy (s*s* +1 + s]s] +1 ) + rs z s z +1 } + n z 
j 

with J x y = Ji + 2J 4 and j* = J l /2 + 5J 4 . 

The effective Zeeman term Tiz and the effective XXZ anisotropy are given 
respectively by 

T _ ~2 



H 



(7) 
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and 



J z Ji/2 + 5J 4 



A eff = Z- = . (8) 

J*w Ji + 2J 4 V ; 

Thus the problem of magnetized states with m = 1/2 of the original system is mapped 
onto that of nonmagnetic states of the Hamiltonian (|6]). It is well-known for the S — 1/2 
XX Z chain that the ground state has the Neel order for A c fr > 1 and magnetic/non- 
magnetic excitations are gapped. The ordered phase is separated from the gapless one 
by a critical point A e fr = 1; the trasition at the critical point is of the Kosterlitz-Thouless 
type with rj = rf = 1, where rj and rf are critical indices defined by 

(S%S?) ~ (-l)V-" and (S'SZ) ~ (-l)V-^ (r -> 00) (9) 

respectively. Turning back to the original system, this implies that the original ladder 
(H) has a magnetization plateau at m — 1/2 under the condition 



J 4 > - Ji. 
o 



(10) 



The plateau state is Neel ordered in the language of the pseudo-spin S, as shown in 
figure 2. 



/Or 



- J O LJ O L - J O L - J / 



i0\ 



Figure 2. Physical picture of the m = 1/2 plateau state. Open ellipses represent | jf) 
and shadowed ellipses | JJ.) in the language of the pseudo-spin S. 



4. Numerical approach 

The phase boudnary J 4 = (1/6) Jj obtained in the previous section is valid for Ji, J4 <C 1. 
To obtain a phase diagram for wider range of the parameters, we have performed the 
numerical diagonalization of the original Hamiltonian (|I|) by the Lanczos method. 

The most powerful method at the present stage to determine the KT phase 
boundary is the level spectroscopy (LS), which was developed by one of the present 
authors (KO) and Nomura |38|-p|. In the LS analysis the critical point is determined 



from the level cross between the two relevant low-lying excitation gaps with a common 
scaling dimension, that is, the same dominant size correction including the logarithmic 
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one. In the present case, we use the following three gaps 

A 1 = ^[E n (L,M + l) + E 7r (L,M-l)-2E (L,M)] (11) 

K = E 7T (L,M)-E (L J M) (12) 
A^ = E^(L,M)-E (L,M) (13) 

(2) 

where Ek(L,M) and E\. (L,M) respectively denote the lowest- and the second lowest 
eigenvalues of the Hamiltonian Ho with the system size L, (Y^ji^ij + ^ij)) = M. These 
gaps Ai, A,,- and A^ govern the long-distance behaviors the S x — S x , S z — S z and dimer 
correlations, respectively, in the language of the pseudo-spin S. The critical indices rj, 
rf and rf can be connected with Ai and A„- by 

LA, LA. ^ = LAg) 

7TV S TIV S TIV S 

respectively, where v s is the spin-wave velocity estimated by 

v s = (E h (L, M) - E (L, M)) (L -> oo) (15) 
with k\ = 27r/L. 

The behaviors of r] and vf near the KT critical point for large L are 

»7 = l-^/o 7^ = l-iy (l + 2*) (16) 

respectively, where i denotes the distance from the KT critical point, and y represents 
the lowest order finite-size correction, which leads to the logarithmic size correction of 
the order of 1/ log L at the KT critical point. Thus the KT critical point can be obtained 
from rj = rf within the lowest order of the finite-size corrections. As shown in equation 
(|T6|), the logarithmic corrections to rj and rf are serious for finite systems. Then the 
relation 77 = 1 (or rf = 1) is not a good indicator for the KT critical point. 

Figure 3 shows the behaviors of rj and rf as functions of J4 for two cases (a) J\ = 0.5 
and (b) J\ = 1.0 for the length of the ladder L = 16. From the crossing points between 
77 and r/ z , we see that the critical point J 4c = 0.10 for J\ = 0.5, and J 4c = 0.22 for 
J\ = 1.0. The quantity 1 — rj = 1 — rj z at J 4c represents the magnitude of the logarithmic 
size correction yo/2 around L = 8 ~ 16. If we determine the KT critical points from the 
condition 77 = 1, we obtain J 4c = 0.18 for J\ = 0.5, and J 4c = 0.36 for J\ = 1.0, which 
are fairly larger than those by use of the LS. 

To check the KT universality class, we also calculated the central charge c by use 

of 

where e g is the ground-state energy per one rung at m = 1/2 for an infinite system. As 
can be seen from figure ^|, the central charge is c = 1 near the KT critical point and 
rapidly decreases in the gapped region, as is expected from |52|,f43|]. This fact supports 
the KT nature of this transition. 
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Figure 3. Critical exponents n and rf as functions of J4 for (a) J% — 0.5 and (b) 
Ji = 1.0 cases. The central charge c (see equation (|l7|)) and the quantity x (see 
equation (|T§|)) are also shown. 



We also checked the KT nature by calculating 

^=^l±pll (18) 

which should be 1/2 in the lowest order of the renormalization equation in the gapless 
region, as was shown by Kitazawa and Nomura [H3 . We can clearly read from figure 
[3] that the quantity x is very close to 1/2 in the gapless region, also verifying the KT 
nature. 

In figure § we show the KT critical points estimated by the LS (i.e. rj = r] z ), as 
well as those obtained by assuming 77 = 1 (conventional method), for the system sizes 
L = 8, 12, 16. We see that the LS result is consistent with that by strong coupling 
approach in §2, while the result obtained by the relation 77 = 1 is not. 



5. Discussion 



In the previous section, we have obtained the phase diagram from the numerical 
diagonalization data through the LS analysis, and shown the KT nature of the transition 
between the plateau and gapless states. Our numerical results are quite consistent with 
those of the strong coupling approach in §3. 

In our previous work |3"2| , where we restricted ourselves to the case of Ji = 1, the KT 
transition was also suggested by the the Roomany-Wyld approximation for the Callen- 
Symanzik /^-function |[45|| . However, the critical exponent rj calculated numerically on 
the basis of the conformal field theory indicated a different value rj ~ 0.5 around 
J 4c ~ 0.05 estimated by the phenomenological renormalization. Since, as already 
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Figure 4. Phase diagram in the magnetized ground state at m = 1/2. Circles and 
diamonds are the critical points determined by LS for L = 8 and L = 12, respectively. 
They are almost independent of the system size. The marks + and x indicate the 
points where 7/ = 1 for L — 12 and L = 16, respectively. They significantly depend on 
L. The result from the strong rung coupling approach J^ c ~ Ji/6 is also displayed. 



shown in $3, the critical exponent rj should be unity at the KT transition point, our 
previous result was not self-consistent. This inconsistency has been completely solved 
by the present analysis. Namely, our previous result J 4c = 0.05 ± 0.04 for J\ — 1 
significantly deviates from the reliable result obtained above. It is consistent with 
the criticism p8| [39| , |4TH that the phenomenological renormalization leads to serious 
underestimation of the gapless region for the KT transition. 

The phase diagram suggests that, for larger J 1; the critical value J 4c becomes larger 
than the limit of a strong rung coupling. It implies that the quantum fluctuation due to 
the leg coupling (Ji) suppresses the plateau formation. In the special case J\ = 1, the 
critical value is J 4c ~ 0.24. Since Sr 14 Cu2 4 04i is expected to have J\ — 1 and J 4 ~ 0.14 
as explained in §1, it is not enough for the appearance of the plateau at m = 1/2. 

An example of the magnetization curve was already shown in our previous work |32| , 
to which readers should refer. At finite temperatures the magnetization plateau will be 
smeared, which remains for future investigations. 

In conclusion, we physically explained the plateau formation scenario of the S = 1/2 
spin ladder with the ring exchange interaction at each plaquette, and also obtained the 
plateau-gapless phase diagram by both analytical- and numerical methods. The present 
phase diagram is more reliable and covers a wider range of parameters than our previous 
one 0. 
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